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Motivated by the recent experimental evidence of commensurate surface 
CDW in Pb/Ge(lll) and Sn/Ge(lll) ^/3-adlayer structures, as well as by the 
insulating states found on K/Si(lll):B and SiC(OOOl), we have investigated 
the role of electron-electron interactions, and also of electron-phonon coupling, 
on the narrow surface state band originating from the dangling bond orbitals 
of the adsorbate. We model the problem by an extended two-dimensional 
Hubbard model at half-filling on a triangular lattice. We include an on-site 
Hubbard repulsion U and a nearest-neighbor V, plus a long-ranged Coulomb 
tail. The electron-phonon interaction is treated in the deformation potential 
approximation. We have explored the phase diagram of the model including 
the possibility of commensurate 3x3 phases, using mainly the Hartree-Fock 
approximation. For U larger than the bandwidth we find magnetic insulators, 
possibly corresponding to the situation in SiC and in K/Si. For smaller U, 
the inter-site repulsion V can stabilize metallic CDW phases, reminiscent of 
the 3x3 structures of Sn/Ge, and possibly of Pb/Ge. 
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I. INTRODUCTION 



Pb and Sn -\/3-adlayer structures on the (111) surface of Ge have recently revealed a 
charge density wave (CDW) instability to a low temperature reconstructed 3x3 phase. 
|l],|2|.|3],[|,[| The low temperature phase is either metallic - as seems to be the case for 
Sn/Ge(lll) - or weakly gapped, or pseudo-gapped, as suggested for Pb/Ge(lll). Related 
systems, like the \/3-adlayer of Si on the (0001) surface of SiC and on K/Si(lll):B, J7J 
show, instead, a clear insulating behavior - with a large gap - with no structural anomalies 
or CDWs investigated so far. Two-dimensional Fermi surface (FS) nesting in the half-filled 
surface states || has been invoked as the driving mechanism for the CDW instability in the 
case of Pb/Ge, Q but excluded for the case of Sn/Ge. |§ 

As all these systems appear to belong together in the same class, we should begin with 
a discussion that in principle encompasses all of them. We believe the following points to 
be of general validity: i) Poor nesting. An unbiased interpretation of the photoemission 
experiments, @,f|,|5] and a close examination of the existing local density approximation 
(LDA) calculations °f the surface ("adatom dangling bond") half-filled band do not 

indicate good nesting of the FS at the surface Brillouin zone (BZ) corner K = (47r/3a, 0). 
We believe this to be equally true for Pb/Ge as for Sn/Ge, contrary to what stated in 
the literature. [|1] Fig. 1, showing the LDA surface band dispersion for the test-case of 
Si(lll) /Si, H as well as the corresponding FS and Lindhard density response function, 
provides a concrete illustration of these statements, ii) Importance of electron-electron 
interactions. The width W of the surface band is relatively small (W ~ 0.5 eV for Pb- 
Sn/Ge, « 0.3 eV for SiC). Moreover, this band is half-filled. These facts call for a careful 
consideration of electron-electron interactions, and not just electron-phonon (e-ph), as a 
possible source of instability. This is reinforced by noting the different phenomenology of 
SiC and K/Si:B with respect to Pb-Sn/Ge, the stronger insulating character of the former 
paralleling closely their stronger electron-electron repulsions, connected with larger bulk 
semiconducting gaps. Hi) Weakness of LDA. LDA is certainly suitable for calculating 
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realistic bands in a weakly interacting system, but is less reliable, at least without including 
spin, in predicting the instabilities of a narrow band. For instance, the phenomenology 
of SiC(OOOl) - suggesting a Mott-Hubbard insulator - is unreproducible by LDA. The very 
onset of a CDW on Sn/Ge(lll) does not seem to be predicted by recent LDA calculations. || 
iv) Different mechanisms for 3x3 CDW instabilities. There are at least two different 
mechanisms which can influence the CDW formation: a) on-site, and nearest-neighbor (n.n.) 
electron-electron repulsion b) on-site effective attraction (negative Hubbard-?/ term) of e-ph 
origin. Of these, the n.n. repulsion naturally suggests, as we shall see, the 3x3 surface 
periodicity, which is found experimentally. Electron-phonon alone would not in itself appear 
to drive a 3 x 3 CDW. At weak coupling, the band structure dominates and incommensurate 
structures could be preferred. At strong coupling, the frustration associated to the triangular 
lattice, will favor, in general, a superconducting ground state over a CDW phase. Jl^ . 



The approach we take here is based on an extended Hubbard-Holstein model. It is by 
necessity a "non-first-principle" approach and, as such, it has no strong predictive power. 
We have found it very helpful, however, in clarifying the possible scenarios as a function of 
physical parameters. 

II. MODEL 

The basic starting point is the half-filled surface state band which one calculates in 
LDA, and is found to lie in the bulk gap. [l|,|J We write our effective Hamiltonian as follows: 

BZ 



k cr 



H = X) X) £ k C k a°k,a + H ph + H c _ ph + H mt , (1) 



where a is the Bloch surface state, whose Wannier function is centered on the adatom, as 
can be seen in Ref. @. H [nt includes correlation effects which are not correctly accounted 
for within LDA, which we parametrize as follows: 

H int = U^2n r ^n r ,l + ~ K-r'Or ~ 1)(tV ~ 1) • (2) 
r 1 r ^r' 
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Here U is an effective repulsion (Hubbard-^7) for two electrons on the same adatom Wannier 
orbital, and V T - r i is the direct Coulomb interaction between different sites r and r'. We 
have considered two models for V r - r r. 1) a truncation to a n.n. V, and 2) a model with a 
Coulomb tail V r - r > = Va/\r — r'|. The results for case 2) are qualitatively similar to those 
of 1), and will be discussed elsewhere. |TT| LDA estimates of the bare repulsions U and 



V between two electrons on the same and on neighboring Wannier orbitals are - for our 
test case of Si(lll) /Si - of about 3.6 eV and 1.8 eV respectively. || Screening effects by 
the the underlying bulk are expected to reduce very substantially these repulsive energies. 
A conservative lower bound for U and V is obtained dividing their bare values by the 
image-charge screening factor, (e + l)/2 ~ 6, yielding U = 0.6 eV, and V = 0.3 eV. As 
for the e-ph interaction, in principle both the on-site Wannier state energy and the hopping 
matrix elements between first neighbors depend on the positions of the adatoms. Within the 
deformation potential approximation, we consider only a linear dependence of the on-site 
energy from a single ionic coordinate (for instance, the height z r of the adatom measured 
from the equilibrium position), and take H e _ p ^ = —gJ2 r z r( n r ~ !)■ The free-phonon term 
will have the usual form = J2k ^Mk (^j^k + 1/2 J, where is the phonon annihilation 
operator. An order-of-magnitude estimate for g is ~ 1 eV/A. 

III. PHASE DIAGRAM. 

We consider first the purely electronic problem in the absence of e-ph interaction. We 
start the discussion from particular limiting cases for which exact statements, or at least 
well-controlled ones, can be made. 

Large positive U : the Mott insulator. For U 3> V, W, the system is deep inside the 



Mott insulating regime. |12| Within the large manifold of spin degenerate states with exactly 
one electron per site, the kinetic energy generates, in second order perturbation theory, the 
Heisenberg spin- 1/2 antiferromagnet as the effective Hamiltonian governing the spin degrees 
of freedom, H cS = Y^,(ij) JijS ri • S rj , with = 4\tij\ 2 /U. |T2[ For our test case of Si(lll)/Si, 
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the values of the hoppings are such that J2/J1 ~ 0.12 while the remaining couplings J3, • • • 
are very small. Antiferromagnetism is frustrated on the triangular lattice. Zero temperature 
long range order (LRO) of the three-sublattice 120°-Neel type - a commensurate spiral spin 
density wave (s-SDW) - is nevertheless likely to occur for the spin 1/2 antiferromagnetic 
Heisenberg model with n.n. coupling J\, and for small enough J 2 . In summary, we expect 
for large values of U a wide-gap Mott insulator with a s-SDW (spins lying in a plane, forming 
120° angles), and a 3 x 3 magnetic unit cell. This is, most likely, the state to be found on 
the Si-terminated SiC(OOOl) surface at T=0. 

Strong inter-site repulsion: an asymmetric CDW with three inequivalent 
sites. The e-ph coupling can effectively reduce U, but not V. Therefore, it is of interest to 
consider the hypothetical regime W < U <C V. In order to minimize the interaction energy, 
the system will prefer a 3 x 3 CDW with two electrons on one sublattice (A), a single electron 
on another sublattice (B), and zero electrons on the third sublattice (C). (See Fig. 2.) The 
spin degeneracy associated with the (unpaired) electron on sublattice B can be removed, 
owing to £2, which leads to an effective spin- 1/2 Heisenberg Hamiltonian within sublattice B, 
with a weak antiferromagnetic exchange constant J = At\/U. JX1|J . Summarizing, we expect 
in this regime a strong 3x3 asymmetric CDW with three inequivalent sites (a-CDW), and 
a spiral 3\/3 x 3-\/3 SDW, governing the unpaired electron spins, superimposed on it. This 
a-CDW is not compatible with the experimental findings on Pb-Sn/Ge, but it could be 
realized in some other case. 

Mean-field theory. In order to get an idea of any additional phases for smaller U, and 
of the possible phase diagram of the model we turn to a mean field analysis. The first issue 
is the possibility of magnetism. For small values of the interactions U and V, the Stoner 
criterion can be used to study the instability of the paramagnetic metal obtained from 
LDA calculations. The charge and spin susceptibility are given, within the random phase 
approximation, by xc(q) = 2x„(q)/[l + (U + 21/ q )x (q)], and xs(q) = Xo(q)/[l - ^Xo(q)], 
where Xo is the non-interacting susceptibility. The divergence of xs is governed, in this 
approximation, by U only. Since Xo(q) is finite everywhere, a finite U is needed in order 
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to destabilize the paramagnetic metal. The wavevector q* at which \s first diverges, by 
increasing U, is in general incommensurate with the underlying unit cell. The instability 
is towards an incommensurate metallic spiral SDW. |13| As for the charge susceptibility, a 
divergence can be caused only by an attractive Fourier component of the potential V^. V q has 
a minimum at the BZ corners ±K, with V±k — — 3V. (V±k ~ —1.54V if a Coulomb tail is 
added). This minimum leads to an instability towards a 3x3 CDW as (£7+2Vk)x (K) = — 1. 

In general, the small coupling paramagnetic metal is surrounded by an intermediate cou- 
pling region, where complicated incommensurate - generally metallic - solutions occur. For 
stronger U and V, commensurate solutions are privileged. [|n| In view of the fact that a 3 x 3 
CDW is experimentally relevant, we concentrate our analysis on the simplest commensurate 
phases. These are easy to study with a standard Hartree-Fock (HF) mean-field theory. (11 



In particular, we restrict ourselves to order parameters associated with momentum space 
averages of the type (c^^c^y) and (c^ a Ck±j£ )(T >) , i.e., the uniform magnetization density m, 
the K-component of the charge density pk = (1/^0 Ek,u( c k cr c k-K,o-), and the K-component 
of the spin density Sk = (1/2 A/") Z)k Sa,/3( c iL a(^)a/3 c k-K,/3)- Pk and Sk have phase freedom, 
and are generally complex: pk = \pk\^ p , etc. The role of the phase is clarified by look- 
ing at the real-space distribution within the 3x3 unit cell. For the charge, for instance, 
{ih^j) = 1 + 2|pk| cos (27rpj/3 + where pj = 0, 1, 2, respectively, on sublattice A, B, and 
C. The e-ph coupling is included but, after linearization, the displacement order parameter 
is not independent, and is given by (zjg) = (g/Mu^PK- Only the phonon modes at ±K 
couple directly to the CDW. The phonon part of the Hamiltonian can be diagonalized by 
displacing the oscillators at ±K. This gives just an extra term in the electronic HF Hamil- 
tonian of the form A£/(p K pK + H.c), with an energy AU = — g 2 /Mu; K which is the relevant 
coupling parameter. This term acts, effectively, as a negative-f/ contribution acting only on 
the charge part of the electronic Hamiltonian. 

The mean-field solutions must be compatible with the symmetry of the problem. A 
symmetry analysis of the Landau theory f|TJ] built from the order parameters m, pk, and 
Sk shows that: 0] i) A CDW can occur without concomitant magnetism, ii) A SDW can 
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occur without a CDW only in the form of a 120° spiral SDW. In all other SDW 
implies also a CDW. Hi) The simultaneous presence of a SDW and a CDW implies, generally, 
a finite magnetization m, unless the phases of p and S are such that 2<p p + (f) a = tt/2 + mr. 

We present a brief summary of the mean- field HF calculations for arbitrary U, V, and g. 
The main phases present in the HF phase diagram are shown in Fig. 2 for the case of g = 0. 
(The diagram is qualitatively similar for g ^ 0.) The s-SDW dominates the large U, small 
V part of the phase diagram, as expected from the Heisenberg model mapping at U — > oo. 
This is the Mott insulator phase, probably relevant for SiC. Its HF bands are shown in Fig. 
3(a). There is however another solution of the HF equations in the large U, small V region. 
It is an insulating state characterized by a linear z-SDW plus a small CDW with <p p = 0, 
accompanied by a magnetization m z = 1/3. It lies above the s-SDW by only a small energy 
difference (of order t x per site), and it could be stabilized by other factors (e.g., spin-orbit). 
A recent LSDA calculation for v3-Si/Si(lll) has indicated this z-SDW as the ground state, 
at least if spins are forced to be collinear. |J The HF bands for this solution are shown in 
Fig. 3(b). By increasing V, the energies of the s-SDW and of the z-SDW tend to approach, 
until they cross at a critical value V c of V. We find V c /t\ « 3.3 at U/t\ = 10. As V > V c , 
however, the insulating a-CDW prevails. For small values of U and V, or for large enough 
e-ph coupling g, a metallic CDW with (p p = (m-CDW) is found. This may be relevant 
for the case of Pb and of Sn/Ge(lll). (See Fig. 3(c) for the HF bands.) The degree of 
metallicity of this phase is much reduced relative to the undistorted surface (pseudo-gap). 
We note that the e-ph interaction can stabilize the <f) p = m-CDW also at relatively large 
U, by countering U with a large negative AU. With g = I eV/A, M Si = 28, and c^k ~ 30 
meV we get AU ~ — 3t±, sufficient to switch from a s-SDW ground state to a m-CDW for, 
say, U/t 1 = 8 and Vjt\ = 2. Much larger g's would eventually stabilize a superconducting 
ground state. ||TT[ 

In conclusion, we have found within a single phase diagram three phases, the s-SDW, the 
z-SDW, and the m-CDW, which may be relevant, respectively, to SiC(0001), to K/Si(lll):B, 
and to Pb or Sn/Ge(lll). 
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Figure Captions 

Figure 1 Surface state dispersion for Si(lll) /Si, as obtained from LDA (solid squares). 
The solid line is a tight-binding fit obtained by including up to the sixth shell of 
neighbors, t 1: ■ ■ ■ , t 6 . The fit gives h = 0.0551 eV, and t 2 /ti = -0.3494, t 3 /h = 0.1335, 
t 4 /tx = -0.0615, t 5 /h = 0.0042, t 6 /ti = -0.0215. Upper inset: The Fermi surface of 
the half-filled surface band. Notice the quite poor nesting at the BZ corner wavevector 
K = (4-7r/3a, 0). Lower inset: The zero temperature Lindhard function x D (q) f° r the 
half-filled surface band. Notice the two peaks located at qi 0.525K and q 2 1.32K, 
and no feature whatsoever at K. 

Figure 2 Schematic Hartree-Fock phase diagram of the U — V model, at zero electron- 
phonon coupling, for the band structure shown in Fig. 1. Only the most important 
commensurate 3x3 phases have been studied. 

Figure 3 Plot of the HF bands along high symmetry directions of the BZ for the s-SDW 
and two CDW 4> p = solutions: (a) at U/ti — 9 and V/ti = 2, the insulating s- 
SDW (ground state); (b) at U/ti = 9 and V/t\ = 2, the insulating solution with a 
small CDW and m z = 1/3 (meta-stable, the actual ground state being the s-SDW); 
Solid and dashed lines denote up and down bands, respectively, (c) at U/t\ =4 and 
V/ti = 2, the metallic solution with a large CDW and no magnetism. 
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